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摘要 : YORP 效应 由 不 规则 形状 小 行星 在 太阳 辐射 下 的 局 部 反 冲 力 差异 导致 的 净 力 矩 所 引起 ， 
长 期 作用 可 以 改变 小 行星 自转 状态 。 从 小 行星 YORP 效应 的 科学 意义 和 研究 概况 出 发 ， 曾 述 了 
YORP 效应 对 小 行星 演化 的 重要 性 ， 并 详细 介绍 了 小 行星 YORP 效应 的 理论 基础 和 观测 研究 现 
状 。 目 前 仅 有 6 颗 小 行星 的 YORP SRAM TER, 计算 得 到 由 YORP 驱动 的 自转 加 速度 ， 
均 处 于 加 速 自 转 状 态 ， 未 来 有 15 颗 小 行星 可 能 通过 实测 测 得 YORP 驱动 自转 加 速度 。 最 后 探讨 
了 利用 小 行星 光 变 等 数据 筛选 可 实测 候选 体 的 方法 。 
X dé dE: 小 行星 ，YORP 效应 ， 形 状 反 演 模型 ， 热 物理 模型 
中 图 分 类 号 : P185.7 文献 标识 码 : A 


小 行星 蕴含 了 太阳 系 形成 早期 的 信息 ， 是 研究 行星 演化 的 活化 石 。 在 小 行星 演化 史 中 ， 
YORP 效应 扮演 了 重要 角色 。 小 行星 的 自转 特性 统计 研究 表明 ”， 直 径 小 于 40 km 的 小 行星 
的 自转 速率 呈现 非 麦 克 斯 韦 分 布 ， 尤 其 是 近 地 小 行星 ， 呈 现 出 明显 的 双 峰 结构 ， 这 与 基于 碰 
撞 的 模型 得 出 的 分 布 不 一 致 ，YORP 效应 为 这 种 自转 分 布 提供 了 一 种 合理 的 解释 。2011 年 ， 
= Scheeres 等 人 ”提出 YORP 也 是 改变 小 行星 大 小 分 布 的 一 种 有 效 机 制 ， 2014 Æ, Jacobson 
o 等 人 的 研究 表明 ， 目 前 主 带 小 行星 的 大 小 分 布 无 法 以 单纯 碰撞 模型 解释 ， 而 在 模型 中 加 入 
YORP 效应 后 与 实测 结果 非常 吻合 。 因 此 ， 深 入 研究 YORP 效应 是 掌握 小 行星 演化 规律 的 
重要 途径 ， 也 为 理解 太阳 系 演 化 提供 了 重要 依据 。 

近年 来 ， 小 行星 的 空间 探测 使 YORP 效应 的 实测 成 为 可 能 ， 抵 近 探 测 提供 了 小 行星 
的 大 小 、 形 状 、 自 转 等 高 精度 测量 数据 ， 极 大 地 促进 了 小 行星 YORP 效应 的 研究 。2003 


收 稿 日 期 : 2020-11-04; 修 回 日 期 : 2021-02-07 


ARMA: 中 科 院 先导 B(XDB41000000); 国家 自然 科学 基金 (11633009); 空间 碎片 与 近 地 小 行星 防御 科研 项 
(KJSP2020020204, KJSP2020020102); 民用 航天 预 研 项 目 (D020304, D020302) 和 小 行星 基金 会 资助 


通讯 作者 : Wi, meteorzh@pmo.ac.cn 


ChinaXiv 合 作 期 刊 


278 天 文学 进展 39 卷 


^E, 日 本 的 华 鸟 号 (Hayabusa II) 在 对 近 地 小 行星 (25143) Itokawa 探测 过 程 中 ， 成 功 地 
测 得 了 Itokawa 受 YORP 驱动 的 自转 加 速度 ”"。2019 年 ， 美 国 的 冥王 号 (origins spectral 
nterpretation resource identification security regolith explorer, OSIRIS-REx) 计划 也 直接 测 
量 了 (101955) Bennu 受 YORP 驱动 的 自转 加 速度 ”。 
众所周知 ， 近 地 小 行星 对 地 球 和 人 类 生存 环境 构成 威胁 ， 而 近 地 小 行星 的 起 源 演化 和 高 
精度 撞击 风险 也 与 YORP 效应 密切 相关 。Yarkovsky 效应 和 YORP 效应 (包括 轨道 YORP 
效应 ) 是 小 尺寸 的 小 行星 从 小 行星 主 带 迁移 到 近 地 轨 道 的 重要 机 制 。 在 研究 (99942) Apophis 
碰撞 地 球 风险 时 ，Vokrouhlicky 等 人 " 发现 Yarkovsky/YORP 驱动 下 短 时 标的 小 行星 轨道 
改变 量 已 经 可 以 通过 目前 的 天 体 测 量 手段 直接 测量 。 因 此 ， 精 确 测量 潜在 威胁 小 行星 的 轨道 
演化 时 ，Yarkovsky/YORP 效应 已 经 不 可 忽略 。 
YORP 效应 研究 可 以 对 小 行星 的 密度 、 内 聚 力 等 物理 特性 进行 约束 。2018 Æ, Scheeres 
等 人 "研究 了 在 YORP 效应 作用 下 碎 石 堆 结构 (rubble pile) 小 行星 瓦解 与 小 行星 密度 、 内 
聚 力 等 物理 特性 的 相关 性 ， 并 给 出 了 不 同 物理 特性 相关 的 瓦解 时 标 。 研 究 小 行星 双星 系统 的 
YORP 效应 ， 可 以 约束 碎 石 堆 小 行星 的 结构 强度 。 通 过 YORP 效应 来 约束 小 行星 物理 参数 
的 研究 将 得 到 越 来 越 多 的 应 用 。 
YORP 效应 也 可 能 是 小 行星 表面 物质 迁移 的 机 制 之 一 。2018 年 ，Kevin 等 人 "认为 
YORP 效应 是 小 行星 表面 结构 重 塑 的 有 效 机 制 ， 并 模拟 了 小 行星 光谱 斜率 (slope) 的 演化 历 
TE RT Q 型 小 行星 丰 度 高 的 原因 。 在 YORP 效应 作用 下 小 行星 加 速 自转 ， 中 纬度 区 
的 物质 向 交道 区 迁移 并 形成 “赤道 将 ”(equatorial ridge)， 例 如 Bennu 和 Ryugu。 2020 年 
Cheng 等 人 "通过 动力 学 模拟 指出 YORP 效应 是 小 行星 “赤道 俏 ” 形 成 的 关键 机 制 。 但 是 ， 
同年 Michel 4& A" RI Bennu 的 赤道 区 存在 古老 的 陨石 坑 ， 表 明 “赤道 产 ” 可 能 早 在 陨石 
坑 形 成 前 就 存在 了 。 综 上 所 述 ，YORP 效应 在 小 行星 表面 形状 的 演化 中 到 底 扮演 着 什么 样 
的 角色 ， 仍 需要 进一步 研究 。 
YORP 效应 是 活动 小 行星 (active asteroid) 的 活动 性 触发 机 制 之 一 。2019 年 ， 小 行 
Æ (6478) Gault 突然 出 现 两 条 尾 迹 ，Kleyna 等 人 ”通过 尘埃 动力 学 模型 (finson-probstein 
pproach) 分 析 了 该 小 行星 两 条 尾 迹 的 粒 径 大 小 分 布 ， 指 出 活动 小 行星 Gault 的 人 尘埃 喷射 是 
H YORP 效应 导致 的 。2019 年 3 月 ， 冥 王 号 近 上 距离 探测 到 了 Bennu 小 行星 的 表面 粒子 喷射 
PT 事件 ， 将 Bennu 归 为 活动 小 行星 ， 目 前 正在 深入 研究 YORP AUS filio Y 09:593 OAL 。 
2019 年 ，Veras 等 人 提出 YORP 效应 可 能 在 恒星 -行星 系统 的 演化 中 起 到 关键 性 作用 ， 
特别 是 在 主 序 星 至 白矮星 演化 阶段 中 必 经 的 巨星 文 阶段 ，YORP 效应 使 得 较 小 的 行星 和 小 
行星 瓦解 成 “YORP 碎片 盘 ”( 距 恒星 中 心 2~100 AU) 宇 。 这 种 效应 在 巨星 支 恒星 的 小 行星 
盘 的 演化 中 是 普遍 现象 。2020 年 ，Veras 5E A" 研究 认为 白矮星 的 金属 污染 也 是 由 小 行星 分 
裂 引 起 的 ， 并 计算 了 YORP 盘存 在 的 时 标 ”。 关 于 YORP 效应 在 恒星 -行星 系统 的 演化 中 
扮演 的 角色 ， 目 前 仅 有 少量 针对 巨星 支 恒星 外 侧 小 行星 的 研究 ， 未 来 将 成 为 恒星 系统 演化 的 
又 一 个 重要 方向 。 
近年 来 YORP 效应 结合 小 行星 观测 数据 的 研究 有 较 大 进展 ， 本 文 将 对 其 重点 介绍 。 本 
文 第 2 章 主 要 介绍 YORP 效应 的 理论 基础 ; 第 3 章 结合 小 行星 实例 介绍 YORP 效应 的 直 
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接 探 测 ， 第 4 章 介 绍 YORP 效应 下 的 小 行星 自转 周期 和 自转 轴 指 向 的 统计 分 布 研究 ， 并 介 
一 种 参数 化 的 统计 方法 ; 第 5 章 讨 论 了 YORP 筛选 模型 以 及 小 行星 表面 微 结构 、 热 红外 
光束 效应 (thermal-infrared beaming effect) 和 全 球 自 加 热效应 (global self-heating effect) 对 
YORP 效应 的 影响 ; 最 后 对 全 文 进行 了 总 结 和 展望 。 


MR 


2 YORP 效应 的 理论 基础 


2.1 YORP 效应 的 提出 
1901 年 ， 俄 罗斯 工程 师 Yarkovsky 发 现 加 热 自转 的 小 行星 可 产生 沿 迹 方向 的 加 速度 ， 
能 抵消 “以 太 假说 ”中 的 阻力 ”， 即 小 行星 吸收 太阳 辐射 后 表面 升温 ， 继 发 的 热 辐 射 产生 
Wis ! 力 ， 从 而 改变 小 行星 轨道 。 该 效应 被 命名 为 Yarkovsky UM. 1954 E, Radziecskii 
等 人 ™ 指出 辐射 压 是 改变 小 行星 自转 速率 的 一 种 机 制 。1969 年 ，Paddack 发 现 太阳 辐射 对 
MTB MY Ee AGMA, Be “LAREN” (windmill effect) Ree MT BAY 
转速 率 ”。 O'Keefe 和 Paddack 用 流水 冲击 不 规则 形状 石头 的 实验 验证 了 该 效应 。 为 了 纪 
& Yarkovsky, O'Keefe, Radzievskii 和 Paddack 这 四 位 科学 家 在 这 种 非 引力 效应 研究 中 的 贡 
献 ， 该 效应 被 命名 为 “YORP 效应 ” 由 于 对 小 行星 物理 特性 的 认识 不 足 ， 很 长 一 段 时 间 内 
YORP 效应 研究 处 于 停滞 状态 。 直 到 2000 年 ，Rubincam 等 人 ™” 才 首次 对 小 行星 的 YORP 
效应 进行 了 定量 研究 ， 并 给 出 了 YORP 效应 力矩 的 计算 公式 ， 从 而 完成 了 YORP 效应 的 理 
论 体系 的 建立 。2007 年 ，Lowry 等 人 开展 了 小 行星 (54509) 2000 PHS 的 自转 特性 的 研究 
发 现在 YORP 效应 作用 下 其 自转 速率 发 生 了 明显 的 变化 ， 首 次 通过 光 变 直接 测量 了 YORP 
旋转 加 速度 。 这 是 YORP 研究 史上 的 重要 里 程 碑 ， 因 此 小 行星 (54509) 2000 PH5 被 命名 为 
“YORP” 从 此 小 行星 的 YORP 效应 研究 进入 了 实测 时 代 。 
2.2 YORP 效应 的 原理 和 假设 
YORP 效应 原理 如 图 1 所 示 ， 小 行星 的 不 规则 形状 导致 其 对 太阳 光 的 反射 和 热 辐射 的 
再 发 射 局 部 不 均匀 ， 这 种 不 均匀 的 反射 和 热 辐射 再 发 射 的 反 冲 力 在 小 行星 上 施加 了 一 个 净 
力矩 ， 虽 然 该 净 力 矩 很 小 ， 但 长 期 累积 作用 可 使 小 行星 的 自转 状态 发 生 明 显 的 改变 。 小 行星 
E Et geo BAR s FENE, IEEE GAA RAI ARTE FREU UST 。 
= 一 般 来 说 ， 研 究 小 行星 的 YORP 效应 的 过 程 中 ， 基 于 两 个 重要 的 基本 假设 ， 一 是 
Rubincam 提出 的 零 热 弛 耶 近 似 ， 二 是 小 行星 绕 惯 量 主轴 旋转 宇 。 零 热 弛 耶 近 似 对 于 米 级 及 
以 下 的 小 行星 是 不 适用 的 ， 因 为 热 惯 量 的 影响 较为 明显 。 小 行星 绕 惯 量 主轴 旋转 假设 则 是 普 
遍 适 用 的 。 对 处 于 非 惯 量 主 轴 旋 转 状 态 和 轨 旋 共振 下 的 小 行星 ， 其 YORP 效应 不 再 是 一 个 
重要 的 因素 ， 不 做 讨论 ”。2.3 节 将 详细 介绍 不 规则 形状 小 行星 的 YORP 力矩 计算 模型 。 
2.8 YORP 效应 的 力矩 计算 
YORP 效应 对 小 行星 的 形状 十 分 敏感 ， 建 立 精确 的 小 行星 形状 模型 是 计算 YORP HE 
的 前 提 。 通 常用 三 角 齐 分 法 来 构建 小 行星 形状 ， 计 算 各 个 三 角 剖 分 面 元 的 力矩 来 计算 小 行星 
YORP 力矩 。 
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1 ”不 规则 小 行星 表面 的 YORP 原理 示意 图 


Breiter 4& A" 在 局 部 太阳 参考 系 (local solar frame, LSF) 中 给 出 了 YORP 力矩 的 计算 
公式 ， 以 面 元 中 心 为 原点 ， 轴 > 为 面 元 外 法 线 方向 ， 轴 x 指向 子午 线 穿 过 太阳 和 视界 面 的 
交点 方向 ， 轴 Y 则 与 其 组 成 右手 正 交 坐标 系 。 定 义 s 为 指向 太阳 的 单位 向 量 ，m 为 指向 天 
顶 的 单位 向 量 。 对 于 面 元 ?， 入 射 通 量 为 四， 即 太阳 热 辐 射 被 面 元 吸收 的 部 分 。 定 义 辐 照 
J£ (irradiance) 为 一 个 面 元 的 任意 方向 的 入 射 通 量 和 面 元 面积 的 比值 ， 则 面 元 i 的 辐 照 度 
A. 


di 
~ d$ 
考虑 地 形 起 伏 可 能 造成 阴影 ， 引 入 可 视 函 数 (visibility function) v， 若 面 元 之 间 无 遮挡 ， 取 
值 为 1， 反 之 为 0; J 则 是 准 直 辐射 密度 。 在 面 元 上 的 反射 光线 的 向 量 o 表示 为 : 


| i-us] 
o= , 


五 (5) =vJs-n. (1) 


/1— u? sind 


u 


(2) 


其 中 ， 为 天 顶 距 余 强 ，9 为 方位 角 。 而 反射 光 从 面 元 dS 上 以 立体 角 dQ 反射 ， 则 反射 辐 
射 率 ( reflected radiance) L, 为 : 


d$, 
(6) = aaa 3) 
引入 双向 反射 分 布 函数 (bidirectional reflectance distribution function) f,(s,o), "Jf: 


L.(o) = f. (s,o) E(s) . (4) 
对 于 太阳 系 的 热 辐射 的 再 发 射 L:， 同 样 可 表示 为 : 
d?o 
fu es (5) 


X pdSdQ 


chinaXiv:202306.00356v1 


ChinaXiv 合 作 期 


3 期 田 君 ， 等 : 小 行星 YORP 效应 的 观测 研究 现状 281 


W Le 和 Lt 的 方向 一 致 ， 则 可 得 : 


_ L,(o) 
e(o) E 7b(o) E (6) 
Ly(o) 为 黑体 发 射 辐射 率 ， 所 以 Lilo) 表示 为 : 
Li(0) = ort (7) 


可 得 单位 面 元 4S、 单 位 立体 角 a 上 的 辐射 反 冲 力 为 : 
F — d'éA-ó)o p 


daR ando — qu s (8) 
用 R2, 表示 向 阳 面 ，9 表示 整个 表面 ， 对 向 阳 面 积分 ， 单 位 面 元 的 力 的 表达 式 为 : 
dF E 20T^n f! 
I Ze) ], vte moin - = f e(u)ap . (9) 
总 力矩 为 : m 
M= fe s (10) 


其 中 , r 是 由 小 行星 质心 指向 面 元 dS 的 位 置 矢量 。 由 式 (10) 可 知 ， 小 行星 的 力矩 主要 受 代 
表 形 状 的 r+ 和 与 辐射 流量 、 小 行星 热 物理 特性 相关 的 和 两 个 因素 影响 。 小 行星 的 形状 有 较 
RONRIREBUS T, JESUS HER 质点 群 法 等 和 ， 因 此 r 的 计算 容易 ， 而 与 小 行星 热 
物理 特性 相关 的 — 计算 误差 较 大 。 
YORP 效应 影响 小 行星 自转 速率 、 自 转轴 倾角 和 进 动 ， 其 中 对 进 动 的 影响 小 ， 通 常 不 
考虑 。2007 年 ，Nesvorny 等 人 ”给 出 YORP 效应 改变 小 行星 自转 速率 和 自转 轴 倾 角 的 计 
算 公 式 ， 建 立 以 小 行星 质量 中 心 (center of mass, COM) 为 原点 的 本 固 坐 标 系 ， 轴 z 为 其 旋 
转 主轴 ， 轴 z 与 最 小 转动 惯量 轴 一 致 。r 和 %n 定义 不 变 ， 将 指向 太阳 的 向 量变 为 wo，0 All o 
分 别 为 小 行星 的 经 度 与 纬度 。 对 于 忽略 热 惯 量 的 小 行星 ， 总 力矩 为 : 


M= -a | ase xn)(n-no), (11) 
S 
其 中 ， Q = 2Fo(1 一 Dv)/(3ve)， Fo 是 小 行星 处 的 太阳 辐射 流量 ， pv 是 反照 率 ， Vc 为 光速 ， 
不 考虑 表面 谈 挡 影响 。 令 nds — N06 — Ntdoysinb)， 其 中 切 向 向 量 N = u " x 则 式 
(11) 改写 为 : 
M = a | Sr N)max(0,n-M9) (12) 
= EFT max(0,n : no) . 


令 r=ro(1+sR)， 其 中 < 尺 用 来 描述 表面 形状 与 标准 球面 的 偏差 ， 则 有 : 


rx N=eror’T , (13) 


C— 
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其 中 ， 
Rọ cos ó cos 0 + Rg sin o sin 0 
T R,sinócos0 — Recosósin0 | ， (14) 
— Rọ sin 0 
Ro 和 Re DHA RE 9/00 上 的 偏 导数 。r x IN. 可 用 泰勒 级 数 展开 为 : 
rx N =er@T(1+2eR) + O(e*) . (15) 
FF} max(0,n .no) = J， 了 表示 面 元 的 辐 照 度 ， 了 表示 面 元 的 平均 辐 照 度 ， 则 : 
I= +e +O(e), (16) 


Ep, DRRR PREE, ch 是 非 球形 状 引入 的 一 阶 修正 项 。 
2007 年 Nesvorny 等 全 重新 梳理 了 小 行星 YORP 力矩 的 计算 公式 ， 其 中 旋转 部 分 的 
YORP JEX: 


: dQ, - dQ i 
Medis i. Solo te iE (Teh + 2RT) ， (17) 


T, 为 下 在 轴 > 方向 的 分 量 。 式 (17) 的 第 一 项 为 一 阶 项 ， 即 ; 


r() 3 df? 
M” = onte f 0 (18) 
AN , Io 独立 o HT, 一 —R, sin, 则 计算 可 得 : 


T 27 
MY? = -arje | dg 人 doo (19) 
所 以 YORP 力矩 的 旋转 分 量 的 一 阶 项 消失 了 ， 同 理 二 阶 项 的 第 二 项 也 为 0， 故 YORP 
力矩 的 旋转 分 量 为 ; 


= = (2 " 
ke Sa 上 UM E. (20) 


YORP 力矩 的 旋转 分 量 是 一 个 二 阶 项 ， 故 YORP ALKA “ZMIE ZU" m. 
小 行星 的 YORP 自转 加 速度 可 写 为 ™” 


dw M; 


Pra E (21) 
YORP 效应 作用 下 的 小 行星 自转 轴 倾 角 的 改变 量 为 ; 

de M, 

dt Cw’ i2) 


其 中 ，C 为 主轴 的 转动 惯量 ，w 为 自转 角速度 ，e 为 自转 轴 倾 角 。 值 得 注意 的 是 ， 在 推导 
YORP 力矩 旋转 分 量 时 ， 出 现 的 勒 让 德 二 阶 项 表示 为 : 


P,(cos€) = 53 cos?c— 1). (23) 
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=e 54.7? 和 125.3° 时 ，YORP 力 盾 旋转 分 量 为 0。 这 两 个 自转 轴 倾 角 和 Koronis 小 行星 
族 Slivan 态 ”中 的 轴 指 向 相近 。 

小 天 体 的 轨道 演化 过 程 中 ，Yarkovsky/YORP 效应 与 碰撞 、 引 力 摄 动 、 喷 流 反 冲 效 
应 等 作用 都 需要 考虑 。 表 1 列 出 了 Yarkovsky/YORP 效应 等 非 引力 效应 与 引力 摄 动 的 对 
比 。 近 20 年 来 ， 随 着 小 天 体 的 高 精度 天 体 测 量 和 光度 测量 的 发 展 ， 我 们 进一步 认识 到 
Yarkovsky/YORP 效应 对 它们 运动 状态 演化 的 重要 性 。 


表 1 非 引力 效应 对 10 cm~10 km 的 小 行星 摄 动 量 级 范围 四 


作用 力 径 向 加 速度 切 向 加 速度 
引力 GMo ~1 一 
Yarkovsky/YORP 效应 1077 ~ 107" 1078 ~ 107”? 
辐射 压 1079 ~ 107"4 
Poynting-Roberton 效应 一 10719 ~ 10715 
太阳 风 、 洛 伦 效 力 、 等 离子 体 阻 力 < 10715 < 10715 


2006 Æ, William 等 人 指出， 直径 小 于 40 km 的 小 行星 中 YORP 效应 对 运动 状态 的 
影响 比 碰撞 更 重要 ， 因 为 这 类 小 行星 YORP 演化 时 标 小 于 碰撞 时 标 ”。 统 计 分 析 表 明 ™ ， 
YORP 效应 对 直径 小 于 15 km 的 小 行星 影响 更 明显 。 直 接 探测 小 行星 的 YORP 效应 应 选择 
直径 15 km 以 下 ” ， 且 无 喷 流 以 及 明显 的 彗星 活动 的 小 行星 。 该 尺寸 范围 的 小 行星 的 碰撞 
时 标 为 10 ~100 Ma. 

除去 喷 流 和 碰撞 外 ， 还 有 四 类 非 引力 效应 ， 其 中 辐射 压 是 太阳 光 等 辐射 粒子 撞击 在 天 体 
表面 上 产生 的 压力 ， 作 用 力矩 为 零 ，; Poynting-Robertson 效应 是 由 太阳 辐射 引起 的 一 个 
与 尘埃 粒子 运动 方向 相反 的 作用 力 ， 主 要 对 直径 为 1 yml mm 的 尘埃 粒子 有 影响 ”;， 太 
阳 风 、 洛 伦 兹 力 、 等 离子 体 阻 力 对 5 ~100 um 的 带电 人 尘埃 粒子 有 明显 作用 号 。 如 表 1 中 所 
示 ， 后 两 种 非 引力 效应 对 直径 10 cm~10 km 的 小 行星 影响 比 Yarkovsky/YORP 效应 弱 3~5 
个 量 级 。 尽 管 与 引力 摄 动 相 比 ，Yarkovsky/YORP 效应 作用 力 微弱 ， 但 其 长 期 作用 不 仅 会 
明显 影响 小 行星 的 自转 状态 ， 而 且 主 导 小 行星 族 的 动力 学 扩散 ”， 是 驱动 小 行星 演化 的 主 
要 非 引力 效应 。 

c 2.4 ”其 他 类 型 的 YORP 效应 

© YORP 效应 也 是 双 小 行星 以 及 多 元 小 行星 系统 形成 的 重要 机 制 之 一 。 双 小 行星 在 演化 
过 程 中 受到 的 YORP 效应 ， 被 称 为 BYORP 效应 。BYORP 效应 最 初 是 由 Ćuk 4 AP" 提 
出 。2010 年 ，McMahon 等 人 四 完善 了 探测 BYORP 效应 的 技术 方法 ， 且 提出 了 BYORP 
效应 的 一 个 分 析 模 型 ， 首 次 应 用 于 小 行星 1999 KW4 的 研究 了 。BYORP 分 析 中 将 双 小 行 
星 看 成 一 个 整体 ， 并 计算 双 小 行星 系统 的 平 运动 偏差 : 


3 /d 
aAM = -4^ (5) a (24) 


Pravec"" f£ BAP (binary asteroid parameters) 中 列 出 了 可 能 存在 平 运动 偏差 的 双 小 行 
星 。BYORP 效应 可 以 改变 双星 系统 的 轨道 构 形 。 次 星 和 主星 质量 比 大 于 0.2， 并 且 处 于 轨 
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旋 同 步 的 双 小 行星 系统 中 ，BYORP 效应 会 使 其 收缩 或 者 扩张 。2014 年 Taylor 等 人 ” 指 
出 ， 只 有 当 BYORP 效应 与 潮汐 作用 相互 抵消 时 迁移 才 会 停止 ， 并 且 不 会 有 明显 的 潮汐 耗 
散 。 双 轨 旋 同步 系统 在 BYORP 效应 作用 下 开始 收缩 会 导致 系统 瓦解 ， 从 而 相互 远离 的 情 
况 ， 目 前 还 没有 详细 的 研究 。 次 星 和 主星 质量 比 小 于 0.2 的 双 小 行星 系统 ， 通 常 稳定 态 为 
个 单 同 步 态 ， 即 次 星 处 于 轨 旋 同步 ， 主 星 则 处 于 加 速 旋转 状态 ， 如 果 次 星 的 BYORP 效应 
为 负 ， 则 系统 处 于 收缩 状态 ， 并 持续 很 长 时 间 ; 如 果 次 星 的 BYORP 效应 为 正 ， 则 系统 扩 
张 ， 与 潮汐 牵引 方向 相同 ， 最 后 失去 同步 性 。 除 了 轨道 效应 外 ，BYORP 效应 还 影响 系统 的 
轨道 平面 指向 ， 即 导致 双星 系统 轨道 极 移 ”。 

2014 年 ，Jacobson 等 人 ”指出 ，BYORP 效应 是 形成 宽 异 步 双星 系统 (wide asyn- 
chronous binary) 的 关键 机 制 之 一 ， 在 观测 的 9 个 宽 异 步 双 星系 统 中 有 6 个 双星 系统 中 的 小 
行星 自转 周期 接近 临界 自转 速率 ， 表 明 YORP 效应 可 能 是 宽 异 步 双 星系 统 形成 的 主要 原因 。 
2011 年 ，Jacobson 等 人 ”提出 在 BYORP 和 潮汐 引力 之 间 存 在 趋 近 平衡 的 观点 ， 对 已 知 
BYORP 偏 移 的 双 小 行星 系统 ， 可 通过 平衡 态 来 确定 快速 旋转 主星 内 的 潮汐 耗 散 ， 从 而 求解 
“ 碎 石 堆 ” 结 构 小 行星 的 力学 强度 等 物理 参数 。2011 Æ, Steinberg 等 人 ”提出 BYORP t 
应 和 YORP 效应 的 强度 正 相 关 ， 并 指出 了 双 小 行星 系统 在 BYORP 作用 下 的 轨道 倾角 的 演 
化 过 程 。 

Rubincam 4 A? 提出 “轨道 YORP 效应 ” (orbital YORP effect) 可 导致 小 天 体 轨道 
偏 移 。 对 于 南北 形状 不 对 称 的 小 行星 ， 受 轨道 偏心 率 影 响 ， 一 个 公转 周期 内 反 冲 力 累 积 
和 用 无 法 抵消 ， 从 而 导致 小 行星 半 长 轴 和 偏心 率 发 生 改 变 。 小 行星 的 热 惯 量 会 减弱 轨道 
YROP 效应 ， 反 照 率 越 高 ， 小 行星 的 红外 辐射 就 越 少 ， 受 到 热 惯量 的 影响 就 越 小 ， 轨 道 
YORP 效应 就 越 大 。Rubincam 等 人 ”指出 Apophis 的 轨道 YORP 效应 会 对 小 行星 轨道 产 
生 与 Yarkovsky 效应 同等 量 级 的 影响 。 轨 道 YORP 效应 对 反照 率 大 、 自 转速 度 慢 ， 且 自转 
轴 与 轨道 面 不 垂直 不 平行 的 小 行星 才能 体现 出 影响 。 小 行星 自转 速率 被 YORP 减速 后 ， 轨 
id YORP 效应 的 影响 会 变 得 明显 ， 所 以 在 研究 近 地 小 行星 的 长 期 撞击 风险 时 必须 要 加 以 
考虑 。 

“ 切 向 YORP 效应 ”是 由 小 行星 表面 亚 米 级 凸 起 结构 引起 的 平行 于 表面 的 反 冲 力作 用 ， 
也 称 为 TYORP 效应 ”” 。 表 面 凸 起 结构 较 多 的 小 行星 上 ，TYORP 效应 量 级 可 与 YORP 效 
应 相当 ”。TYORP 效应 会 加 速 小 行星 自转 ， 被 认为 是 YORP 效应 作用 下 缺失 减速 自转 小 
行星 的 重要 原因 。 


3 YORP 效应 的 直接 探测 


本 章 将 介绍 YORP 效应 的 一 种 研究 方法 一 一 直接 探测 法 。 小 行星 自转 轴 倾 角 的 演化 十 
分 复杂 ， 了 且 自 转轴 倾角 的 改变 很 难 直接 探测 到 ， 所 以 直接 探测 法 主要 计算 小 行星 的 YORP 
旋转 加 速度 。 目 前 直接 探测 小 行星 YORP 自转 加 速度 的 方法 依赖 两 种 模型 ， 分 别 为 光 变 曲 
线形 状 反 演 模 型 和 热 物理 模型 。 


T 
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3.1 ”基于 光 变 曲线 的 形状 反 演 模型 

基于 光 变 曲线 构建 小 行星 形状 模型 ， 主 要 使 用 包含 多 次 回归 中 不 同 相 位 角 的 高 密 
度 时 序 测 光 数据 ， 近 年 来 稀 玻 测 光 数据 作为 补充 数据 也 应 用 于 形状 模型 构建 。Durech 
等 人 组建 并 一 直 在 更 新 小 行星 形状 模型 的 数据 库 (database of asteroid models from 
inversion techniques, DAMIT) 。 

小 行星 形状 基于 三 种 基本 模型 ， 分 别 是 三 轴 椭 球 模型 、 凸 沉 模型 和 非 凸 模型 ”。 为 了 
更 好 地 体现 小 行星 a 使 用 基于 高 斯 球 模型 的 三 角 剖 分 面 元 来 构建 形状 模型 。 选 取 
合适 的 散射 模型 ， 通 常用 Lommel-Seeliger 模型 ， 对 每 个 面 元 的 流量 进行 求 和 ， 得 到 人 造 光 
ae tie P" 
应 用 傅 里 叶 分 析 来 得 到 小 行星 的 会 合 周期 Tynde ， 再 根据 会 合 周期 计算 小 行星 的 恒 
星 周期 Trotational? 


NS 


Jnd = Lotani Le “oo = Fel) ) (25) 

H, (Leo 一 m uu e M. pU 

求解 小 行星 的 自转 轴 指 向 时 ， 方法 可 以 通过 光 变 曲线 振幅 ， 利 用 最 小 二 乘法 求 得 极 
点 坐标 ; 另 一 种 方法 是 通过 光 变 周期 求解 自转 轴 方 向 ， 即 利用 恒星 周期 和 反 演 方法 对 自转 名 
指向 进行 扫描 ， 得 到 最 优 坐标 并 在 误差 范围 内 验证 其 可 靠 性 。 

基于 光 变 曲线 的 形状 反 演 模型 可 得 到 小 行星 的 形状 、 人 恒星 周期 和 极点 黄道 坐标 等 重要 
的 物理 参数 。 该 方法 探测 YORP 效应 的 核心 思想 是 利用 反 演 得 到 的 小 行星 形状 分 别 产生 常 
周期 和 周期 线性 变化 下 的 人 造 光 变 ， 并 与 实测 光 变 进行 最 优 拟 合 ， 若 后 者 的 拟 合 情况 明显 优 
于 前 者 ， 那 么 用 后 者 最 优 拟 合 情 况 下 的 周期 线性 改变 量 来 计算 与 之 对 应 的 YORP 旋转 加 速 
度 四 四 

计算 过 程 有 两 种 不 同 的 处 理 方式 。 其 一 ， 设 定 自 转 周 期 ， 计 算 人 造 光 变 与 实测 光 变 之 间 
的 相位 偏 移 。YORP 效应 造成 的 相位 偏 移 yp Ft KR AA: 


| p= jt — Tp)? + (wy -- (t — Ty) , (26) 
= ELH, oft To 时 刻 的 自转 速率 wo 的 估计 误差 ，T 时 刻 任何 的 自转 速率 e 都 是 一 个 非 零 值 。 
"T 其 二 ， 在 形状 反 演 模 型 中 增加 一 个 YORP BA, ARE Aste 5; Sci 638 ZA AY 

相位 偏 移 ， 直 至 人 造 光 变 与 实测 光 变 拟 合 的 x? 无 明显 差异 ，X2 最 小 对 应 的 YORP 参数 即 


为 YORP 旋转 加 速度 ” 。 

本 文 以 (1620) Geographos 小 行星 为 例 ， 在 其 1969 一 2019 年 间 的 118 条 光 变 曲线 中 ， 
od ra ati 74 KA HAE, f SE T. Geographos 小 行星 的 
形状 模型 ， 图 2 是 Geographos 小 行星 形状 模型 的 三 视图 。 小 行星 自转 轴 指 向 的 坐标 为 
(52.2 + 6°, 一 40.0 土 7°)， 恒 星 周 期 为 P= 5.223 319 + 0.000 002 h(To = 2 440 229.0). 
通过 最 小 二 乘法 拟 合 74 条 实测 光 变 曲线 与 最 佳 常 周 期 模型 的 人 造 光 变 曲线 的 相位 差 ， 
得 到 YORP 旋转 加 速度 为 v = 1.195 x 1078 rad/d? ( 见 图 3). 

Durech 等 人 分 别 得 到 常 周 期 和 添加 YORP 参数 后 的 人 造 光 变 曲 线 。 图 4 是 部 分 人 造 光 
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VE: a) 轴 > 方向 ; 


相位 偏 移 /min 


a) b) c) 
b) fi y 方向 ; c) # a 方向 。 
2 ”小 行星 (1620) Geographos 形状 模型 的 三 视图 
40 
© 
25 
10 
-5 
LJ © 
-20 
40 000 45 000 50 000 55 000 60 000 
观测 日 期 /d 
图 3 ”人造 光 变 曲线 和 实测 光 变 曲线 之 间 的 相位 偏 移 的 拟 合 曲线 


1969/10/7.2 1994/10/30.9 


0.2 


0=87° 0,—84? a=9° 
0.4 0.6 0.8 1.0 0 0.2 0.4 0.6 0.8 1.0 
自转 相位 自转 相位 


图 4 ”人 造 光 变 曲 线 与 实测 数据 的 对 比 
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变 曲线 与 实测 数据 的 对 比 ， 虚 线 为 常 周期 情况 下 的 人 造 光 变 曲线 ， 与 实测 数据 存在 明显 的 相 
位 差 ， 实 线 为 YORP 参数 为 v = 1.15 x 10 rad/d? 时 对 应 的 人 造 光 变 曲线 ， 与 实测 数据 拟 
合 较 好 ”。 小 行星 的 光度 由 视界 角 6， 太阳 方位 角 90。 和 相位 角 a 给 出 。 方 法 一 相位 偏 移 拟 
合 耗 时 长 ， 且 会 引入 一 定 误差 ;方法 二 更 高 效 。 

迄今 为 止 ， 只 在 6 颗 近 地 小 行星 上 直接 探测 到 了 YORP 驱动 的 自转 加 速度 ( 见 表 2), 
均 为 加 速 自转 状态 ， 还 没有 发 现 YORP 效应 作用 下 减速 自转 的 小 行星 。 


表 2 6 个 已 被 直接 探测 到 YORP 旋转 速率 改变 量 的 小 行星 


r 


小 行星 名 (1620) (1862) (3103) (25143) (54509) (101955) 
Geographos Apollo Eger Itokawa YORP Bennu 

半 长 轴 /AU 1.245 1.470 1.405 1.324 1.006 1.259 

偏心 率 0.336 0.560 0.345 0.20 0.230 0.204 

绝对 星 等 /mag 15.6 16.3 15.4 19.2 22.7 20.9 

平均 直径 /km | 2.56 +0.05 1.45+0.22 1.80 土 0.30 0.327 0.001 0.113+0.002 0.265 + 0.0375 

几何 反照 率 0.16 +0.01 0.25 士 0.08 0.38 士 0.13 0.345 士 0.002 0.115+0.004 0.044 士 0.002 

v 恒星 周期 /h 5.223 3364 3.065 448+ 5.710 1564 12.132 38+ 0.202 899 41+ 4.297 46+ 


2 x 1078 3 x 1078 7 x 1078 2 x 1074 10 2 x 107? 


YORP (1.15 +0.15) (5.32: 1.3) (1.10.5) (3.54+0.38) (3.50.3) — (5.86 0.91) 
/rad - d^? x10^5 x10^5 x1075 x1075 x107 x1078 
时 标 /Ma 6.9+0.9 25:05 5242.2 > 0.23 0.58 + 0.06 ~15 
光 变 时 间 跨 度 | 1969 — 2008 1980 — 2007 1987 — 2012 2000 — 2007 2001— 2005 1999 一 2018 
参考 文献 [5g] lar, ez] [53 64) 回 lEs] [5]. 66] 


除 此 之 外 ， 还 有 (1865) Cerberus" 和 (2100) Ra-shalom"" T 等 多 个 小 行星 确定 存在 
YORP 效应 ， 但 目前 仅 给 出 了 YORP 旋转 加 速度 的 估 值 范围 。 其 中 (161989) Cacus 小 行星 
虽然 被 明确 探测 到 存在 YORP 效应 “” ， 且 给 出 了 YORP 旋转 加 速度 为 1.9x10-8 rad/d?, 
但 其 过 于 依赖 1978 年 的 光 变 数据 ， 结 果 的 置信 度 不 高 ， 所 以 本 文 未 将 其 归 入 已 测 得 YORP 
旋转 加 速度 的 小 行星 列表 中 。 

YORP 效应 的 直接 探测 不 仅 需 要 高 密度 光 变 曲线 反 演 形状 ， 而 且 需 要 光 变 曲线 有 较 
长 的 时 间 跨 度 ， 所 以 被 探测 到 数目 较 少 ， 表 3 列 出 了 在 未 来 具备 观测 条 件 且 可 能 计算 出 
YORP 旋转 加 速度 的 小 行星 ， 他 们 是 未 来 需要 关注 的 重点 目标 。 

基于 光 变 曲线 的 形状 反 演 模型 是 探测 小 行星 YORP 旋转 加 速度 的 主要 手段 ， 事 实证 明 
该 方法 非常 可 靠 ”。 不 过 该 模型 对 小 行星 表面 的 微 结构 不 敏感 ， 一 定 程度 上 限制 了 YORP 
旋转 加 速度 的 计算 精度 ， 因 此 结合 热 物理 模型 进行 研究 非常 必要 。 

3.2 ” 热 物理 模型 

早期 针对 不 同 目标 建立 的 热 物 理 模 型 有 STM (standard thermal model), FRM (fast 
rotating model), NEATM (near-earth asteroids thermal model)， 均 是 把 小 行星 作为 无 自转 
的 球体 来 研究 的 简化 热 模型 。 为 了 更 好 地 描述 小 行星 表面 的 温度 分 布 情况 ， 进 一 步 发 展 的 
热 物理 模型 (thermal physical model, TPM) 和 先进 热 物 理 模型 (advanced thermal physical 
model, ATPM) 便 应 运 而 生 。TPM 通 过 红外 实测 数据 得 到 小 行星 的 反照 率 、 比 热 容 和 密度 
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表 3 15 个 在 未 来 的 探测 中 可 能 测 得 YORP 旋转 改变 量 的 小 行星 
小 行星 YORP 半 长 轴 ”绝对 星 等 恒星 周期 
/rad- d? /AU /mag /h 
(1685) Toro" ^ ~ 2.3 x 10-9 1368 1423 10.197 82 + 0.000 03 
(162173) Ryugu = 1.190 19.30 7.631 2 + 0.000 2™ 
(161989) Cacus"" ~ 1.9 x 1078 1.123 17.2 3.755 067 + 0.000 002 
(2100) Ra-Shalom"" (—5.98 ~ 2.03) x 1078 0.832 16.05 19.820 0 + 0.000 3 
(1865) Cerberus" ^ < 0.8 x 1078 1.080 16.84 6.803 286 + 0.000 005.7" 
1998 KY26 7 三 1.233 25.5 0.178 358 3 0.000 007™ 
2000 HB24 四 = 0.815 23.7 0.217 6 土 0.000 6 
2000 YA 四 z 2.386 23.7 0.665 8 + 0.000 1 
2001 SQ3™ = 1.110 21.6 0.062 48 + 0.000 05 
2001 AV43 7 = 1.284 24.6 0.170 1+ 0.000 5 
2006 xv z 1.504 24.2 0.082 978 3 + 0.000 000 3 四 
2007 DD” 2 0.987 25.9 0.074 29 + 0.000 07 ™ 
(951) Gaspra” = 2.210 11.46 7.042 027 + 0.000 0057" 
1996 HW1 = 2.045 15.3 8.762 39 + 0.000 05" 7” 


等 物理 参数 。 但 TPM 假 设 小 行星 是 三 轴 椭 球体 ， 在 计算 YORP 效应 时 存在 局 限 性 。2011 
年 ，Rozitis 等 人 "提出 了 先进 热 物 理 模 型 。 
图 5 展示 了 小 行星 表面 热 辐射 的 主要 类 型 。 使 用 热 物理 模型 计算 小 行星 的 YORP 效应 
时 ， 通 常 不 考虑 多 重 辐射 效应 。 
2014 f£, Lowry 等 人 "应 用 先进 热 物理 模型 计算 了 Itokawa 小 行星 的 YORP 自转 加 速 
度 ， 发 现 与 基于 形状 反 演 模型 得 到 的 YORP 旋转 加 速度 有 较 大 偏差 ( 见 图 6)。 
先进 热 物理 模型 并 没有 考虑 “ 热 红 外 束 流 ” 效 应 和 全 球 自 加 热效应 的 影响 。 不 过 热 物 
模型 在 已 知 小 行星 YORP 驱动 旋转 加 速度 后 可 以 约束 其 密度 、 反 照 率 等 物理 参数 。2018 
年 ， 姜 浩 轩 和 季 江 徽 ” 详 细 介 绍 了 小 行星 热 物 理 模型 在 YORP 效应 的 研究 进展 。 
2015 年 ，Sevetek 等 人 ”给 出 了 考虑 小 行星 巨石 结构 产生 的 谈 挡 效应 和 局 部 自 加 热效应 
时 的 先进 热 物 理 模 型 ， 其 中 面 元 的 反 冲 力 计算 公式 为 : 


4 $ _ , i 
df eou 2 | f ORG r E vdr (27) 
ds (BI nir-r? r-r] 


其 中 , r flr 为 表面 不 同位 置 的 位 置 向 量 ，o 是 Stefan-Boltzmann 常数 ，c 为 光速 ，e 为 热 
KNB, u 为 温度 ，n 为 小 行星 外 表面 的 法 线 方向 癌 量 ，a 为 面 元 法 向 与 位 置 向 量 7 的 夹 角 ， 
D, 表示 小 行星 表面 的 边界 ，T" 表示 小 行星 表面 面 元 。 该 模型 仍 没有 考虑 “ 热 红 外 束 流 ” 交 
应 和 全 球 自 加 热效应 的 影响 。 先 进 热 物理 模型 在 计算 YORP 效应 时 需要 进一步 完善 。 
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太阳 


图 5 ”小 行星 表面 热 辐射 的 主要 类 型 


归 一 化 频率 


0.0077 -5 0 5 X107 


YORP fre $4 JE E /rad-d? 
VE: 虚线 为 光 变 观测 值 。 
6 ”使 用 ATMP 在 不 同 粗糙 度 下 YORP 旋转 加 速度 的 概率 分 布 “ 
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4 YORP 效应 的 统计 分 析 


4.1 ”小 行星 自转 速率 的 分 布 

YORP 效应 对 主 带 小 行星 和 近 地 小 行星 的 影响 程度 不 同 。2008 年 ，Pravec 等 人 四 基于 
观测 数据 ， 绘 制 了 直径 在 3~15 km 间 的 268 颗 主 带 及 “ 越 火 ” 小 行星 (MB/MC) 的 自转 速 
率 分 布 直方 图 ( 见 图 T). MB/MC 小 行星 的 自转 速率 分 布 较为 平坦 ， 其 中 慢 自转 的 小 行星 占 
比较 高 。 对 于 近 地 小 行星 ， 自 转速 率 分 布 则 存在 明显 的 双 峰 特征 ( 见 图 8). 


W 


60 60 
40 40 
T: Tg] 
ES Es 
20 20 
0 
01 2 3.4 56 7 8 91011 12 Ü 123 45 6 7 8 9 10 11 12 
自转 速率 /d+ : 转速 率 /d-! 
YE: 其 中 深 色 是 双 小 行星 ， 灰 色 是 单 小 行星 。 FE: 可 以 看 到 其 存在 明显 的 双 峰 结构 ， 即 快速 
虚线 标识 位 置 为 自转 速率 在 1 ~ 9 d-! 范围 转 和 缓慢 自转 的 近 地 小 行星 较 多 ;这 正 是 
内 样本 的 平均 数 密度 。 YORP 效应 演化 作用 的 证 据 。 
图 7 直径 在 3~15 km 范围 内 的 主 带 小 行星 以 及 图 8 直径 大 于 200 m 的 近 地 小 行星 的 自转 速率 分 
“ 越 火 ” 小 行星 的 自转 速率 分 布 图 a” 


2009 4, Rossi 等 人 ”对 近 地 小 行星 自转 分 布 的 数值 模拟 中 ， 发 现在 不 考虑 YORP 效 

应 的 情况 下 ， 模 拟 结 果 和 观测 存在 明显 偏差 ， 而 加 入 YORP 效应 后 得 到 了 符合 观测 的 结果 。 
由 此 可 知 ，YORP 效应 对 小 行星 的 自转 速度 演化 有 非常 明显 的 影响 。 

= 42 ”小 行星 自转 轴 指 向 的 分 布 

2011 Æ, Hanuš 等 人 ”基于 221 个 小 行星 的 形状 数据 ， 统 计 其 自转 轴 指 向 的 分 布 ， 图 

9 为 直径 小 于 30 km 的 小 行星 的 经 纬度 分 布 。 从 图 中 可 看 到 自转 轴 的 纬度 分 布 呈 明 显 的 双 峰 

特征 ， 而 经 度 分 布 则 没有 明显 特征 。 

为 了 进一步 确定 小 行星 自转 轴 指 向 的 分 布 特点 ，Hanug 等 人 创建 了 一 个 包括 但 不 限 

T YORP 效应 、 碰 撞 等 因素 的 小 行星 自转 演化 模型 ， 模 拟 结 果 与 观测 结果 吻合 ， 呈 现 出 双 

峰 特 征 ( 见 图 10). 

小 行星 自转 轴 指 向 的 聚集 是 YORP 效应 驱动 小 行星 演化 的 有 力 证 据 之 一 ， 这 在 小 行星 

族 中 以 族 成 员 自 转轴 指向 的 聚集 (slivan AS) 来 呈现 。2002 年 ，Slivan 等 人 ”发 现 Koronis 小 


T 


r 
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纬度 /(®) 
a) 


没有 明显 的 特征 。 


图 9 ”直径 小 于 30 km 的 小 行星 的 经 纬度 分 布下 


VE: a) 为 小 行星 极点 纬度 分 布 ， 明 显 呈 现 双 峰 特 征 ， 即 小 行星 自转 轴 指 向 聚集 在 高 纬度 位 置 ，b) 为 经 度 分 布 ， 


15 
10 
a 
x 5 
: : -53 -24 0 24 53 
ARJO) ARSC) 
a) b) 


TE: a) 为 直径 小 于 30 km 的 小 行星 ， 通 过 旋转 模型 模拟 得 到 的 纬度 分 布 ，b) 为 实测 数据 的 纬度 分 布 图 。 


10 ”小 行星 极点 纬度 的 模拟 结果 与 观测 到 的 纬度 分 布 


行星 族 的 9 个 族 成 员 的 自转 轴 指 向 对 齐 的 现象 ， 次 年 Vokrouhlicky 等 人 ”指出 该 对 齐 现 象 


与 族群 成 员 所 受 的 YORP 效应 相关 。 


2003 年 ，Vokrouhlicky 等 人 ”指出 主 带 外 侧 可 能 存在 大 量 处 于 Slivan 态 的 小 行星 。 


2015 4E, Vra&til 55 AP 模拟 了 主 带 内 侧 小 行星 的 演化 ， 指 出 Massalia 小 行星 族 可 能 也 处 于 


Slivan 态 ， 次 年 ，Vrastil 等 人 ”通过 数值 模拟 寻找 小 行星 主 带 中 存在 Slivan 态 的 参数 空间 ， 
发 现 位 于 主 带 外 侧 低 倾角 轨道 上 的 小 行星 很 可 能 处 于 Slivan 态 ， 但 目前 仍 没有 被 观测 证 实 。 


实际 上 ， 主 带 外 侧 的 小 行星 受到 YORP 效应 的 影响 已 经 不 再 明显 ， 轨 旋 共 振 可 能 是 
Slivan 态 的 主要 原因 。 


存在 


2016 年 ，Paolicchi 等 人 ”在 小 行星 族 成 员 的 绝对 星 等 与 轨道 半 长 轴 分 布 图 中 发 现 一 
个 空 腔 区 ， 称 为 “YORP HR" (YORP-eye). *YORP 眼 ” 的 大 小 与 小 行星 族 的 年 龄 相关 ， 


新 生 的 族 不 存在 该 结构 ”。 小 行星 族 的 年 龄 越 大 ,“YORP 眼 ” 对 应 的 星 等 值 越 小 ， 因 此 


“VORP 眼 ” 给 出 了 一 种 独立 估计 小 行星 族 年 龄 的 方法 ”。2019 年 ，Marzari 等 人 "构建 了 
一 个 考虑 Yarkovsky/YORP 效应 和 碰撞 影响 的 小 行星 族 演化 模型 ， 较 好 地 还 原 了 小 行星 族 


的 V 型 图 ， 并 指出 “YORP 眼 ” 是 小 行星 族 在 Yarkovsky/YORP 效应 作用 下 演化 的 结 


FH 
人 小。 
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4.3 ”参数 化 YORP 效应 计算 方法 
参数 化 YORP 效应 计算 方法 是 基于 真实 的 小 行星 物理 参数 和 轨道 参数 ， 将 小 行星 形状 
和 热 惯 量 等 影响 用 YORP 系数 代替 的 方法 。 
2009 ££, Rossi 等 人 给 出 了 YORP 效应 影响 小 行星 自转 速率 的 公式 : 
dw G4 
dt — RU Te = 737 d , (28) 
A, Gi 是 修正 的 太阳 常数 ， 约 为 6.4 x 101? kg - km - $72, a, e, p, D 分 别 是 小 行星 的 轨道 
长 径 、 轨 道 偏心 率 、 体 密度 和 有 效 直径 ，Cy 是 YORP 系数 (或 YORP 因子 )。 
2013 年 ，Rozitis 等 人 "给 出 了 YORP 旋转 加 速度 的 探测 条 件 : 
dw | 8nX Cp 
-一 | > 
at |? TE. 
EH, Cp 是 旋转 精度 系数 ，TeAw 是 光 变 数据 的 时 间 跨 度 ，X 是 在 小 行星 上 探测 到 YORP 
效应 的 置信 水 平 。 小 行星 YORP 效应 的 探测 条 件 和 自转 周期 无 关 。 该 参数 化 的 分 析 方 法 ， 
不 仅 可 以 估算 小 行星 YORP 效应 的 大 小 ， 也 可 以 给 定 光 变 时 间 跨 度 下 探测 到 YORP 效应 的 
概率 ， 或 给 定 探测 概率 下 所 需要 的 光 变 时 间 跨 度 。 该 方法 不 适用 于 快速 自转 (自转 周期 小 于 
2.2 h) 和 缓慢 自转 (自转 周期 大 于 20 h) 的 小 行星 ， 而 且 对 单个 样本 的 分 析 结 果 有 较 大 误差 ， 
但 从 统计 角度 来 看 具有 重要 意义 。 


T 


kx 


(29) 


5.1 ”可 实测 YORP 效应 目标 筛选 

随 着 海量 小 行星 高 密度 光 变 数据 和 稀疏 测 光 数据 的 产生 ， 给 建立 筛选 可 实测 YORP 效 
应 的 小 行星 候选 体 的 方法 提供 了 用 武之 地 。 建 立 一 个 普 适 的 筛选 模型 ， 将 大 大 提高 对 小 行星 
YORP 效应 的 探测 效率 。 我 们 以 可 能 存在 且 可 能 探测 到 YORP 自转 加 速度 的 小 行星 为 基本 
目标 ， 从 形状 反 演 模型 所 需 的 光 变 数据 的 条 件 入 手 ， 建 立 YORP 效应 候选 体 得 选 模型 。 

当然 ， 建 立 完 善 的 筛选 模型 仍然 存在 较 多 的 不 确定 性 。 首 先 ，YORP 效应 的 力矩 与 小 
行星 的 直径 、 密 度 、 热 惯量 等 物理 参数 密切 相关 ， 不 同 小 行星 的 热 物理 参数 和 形状 差异 较 
大 ， 甚 至 小 行星 物理 参数 所 参考 的 标准 都 不 统一 。 其 次 ， 稀 玻 测 光 数 据 是 否 真 的 有 助 于 
YORP 效应 的 探测 仍 不 明确 呈 ， 筛 选 模型 是 否 应 该 考虑 大 量 稀 疏 测 光 数据 的 小 行星 的 存在 。 
还 有 ， 处 于 共振 区 的 小 行星 是 否 应 该 作为 YORP 效应 候选 体 的 备 选 目标 ， 等 等 。 
5.2 ” 微 结 构 对 计算 YORP 效应 的 影响 

基于 光 变 的 形状 反 演 模型 缺乏 对 表面 微 结构 (hn Roe vn Be. IUE BS BU PC) 刻画 的 
能 力 。2009 年 ，Statler EAT 首次 注意 到 YORP 效应 对 小 行星 表面 微 结 构 十 分 敏感 ， 甚 至 
可 以 改变 小 行星 YORP 力矩 的 方向 ， 量 化 微 结构 对 YORP 效应 的 影响 程度 变 得 重要 。2008 
年 ， 了 urech 等 人 指出 Apollo 小 行星 的 表面 微 结构 对 小 行星 的 影响 约 为 其 YORP 效应 的 
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Ree 向 不 超过 595". HEA] 2015 4, Golubov 等 人 
行星 的 微 结 构 对 YORP 效应 的 影响 。 
“ 热 红外 束 流 ” 效 应 和 全 球 自 加 热效应 


| | 


II 


才 尝试 用 数学 方法 来 量化 


行星 粗糙 的 表面 、 较 大 的 凹 坑 、 表 面 巨石 产生 的 阴影 均 会 产生 两 种 效应 ， 分 别 是 “ 热 
流 ” 效 应 和 全 球 自 加 热效应 ”。 图 11 和 图 12 分 别 是 两 种 效应 的 示意 图 。 


PBN TCER 不 规则 小 行星 


mm 


AMT EE Qs 


图 11 RER” AARRE” 图 12 BIS ENnBGRSGEARD 


热 红 外 束 流 ” 效 应 是 指 在 某 相位 角 处 ， 观 测 到 小 行星 表面 红外 辐射 陡 增 的 现象 。 若 
表面 的 粗糙 度 ， 辐 射 方 向 通常 重 直 于 表面 。2011 年 ，Rozitis 等 人 "提出 热 物 理 模型 
虑 “ 热 红外 束 流 ” 效 应 。2012 年 ，Rozitis 等 人 四 发 现 ,“ 热 红外 束 流 ”效应 的 影响 
取决 于 热 惯 量 、 反 照 率 等 ， 而 对 小 行星 表面 粗糙 度 非常 敏感 ， 对 YORP 效应 的 抑制 


着 粗糙 度 的 增加 而 增加 。 Rozitis 等 人 四 通过 与 光滑 表 的 小 行星 的 对 比 研究 ， 发 现 


热 红 外 束 流 ” 效 应 最 高 可 将 YORP 效应 减弱 到 一 半 。 这 表明 ， 如 果 应 用 热 物理 模型 时 不 考 


‘红外 束 流 ” 效 应 ， 计 算得 到 的 YORP 自转 加 速度 可 能 存在 很 大 的 误差 ， 甚 至 是 不 准 


局 部 自 加 热效应 是 指 小 行星 未 被 太阳 照 亮 的 凹面 接受 到 了 其 他 面 元 反射 的 光 ， 寻 致 阴 


影 区 域 被 加 热 并 发 出 热 辐射 的 现象 ， 全 球 自 加 热效应 即 整个 小 行星 上 发 生 的 自 加 热 现象 。 在 


处 理 Itokawa 这 样 的 存在 较 大 阴影 区 的 小 行星 时 ， 必 须 考 虑 全 球 

等 人 "指出 全 球 自 加 热效应 会 降低 YORP 效应 对 小 行星 表面 微 结构 的 敏感 度 。 与 粗 
下 热 红 外 束 流 效应 相似 ， 全 球 自 加 热效应 对 YORP 效应 也 有 抑制 作用 ， 影 响 大 小 约 
为 10%。 
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都 同时 存在 这 两 种 效应 ， 他 们 的 共同 作用 对 YORP 效应 的 影响 是 未 来 研究 的 一 个 方向 。 


6 总 结 和 展望 


本 文 主要 介绍 了 小 行星 YORP 效应 的 观测 研究 进展 ， 从 小 行星 YORP 效应 的 基本 理 


论 、 直 接 探 测 以 及 统计 分 析 三 个 方面 进行 了 论述 。 基 于 光 变 数据 的 形状 反 演 模型 和 热 物 理 


模型 是 探测 YORP 旋转 加 速度 的 有 效 手 段 。 形 状 反 演 模型 需要 时 间 跨 度 较 长 的 高 密度 时 
序 测 光 数据 ， 极 大 地 限制 了 对 小 行星 YORP 效应 的 实测 ， 因 此 目前 仅 有 6 颗 近 地 小 行星 的 
YORP 旋转 加 速度 被 直接 探测 。 热 物理 模型 由 于 没有 考虑 “ 热 红外 束 流 ”效应 和 全 球 自 加 
热效应 ， 降 低 了 探测 YORP 旋转 加 速度 的 精度 ， 不 过 热 物理 模型 能 有 效 地 约束 已 知 YORP 
旋转 加 速度 的 小 行星 的 热 惯量 等 热 物理 参数 。 由 于 小 行星 的 形状 数据 和 热 物理 模型 仍 十 


分 有 限 ， 本 文 还 介绍 了 一 种 参数 化 的 分 析 方 法 ， 将 小 行星 的 形状 和 热 惯量 等 


热 物 理 参数 对 


YORP 效应 的 影响 用 YORP 系数 来 表示 ， 不 仅 有 利于 小 行星 YORP 效应 的 统计 研究 ， 还 可 


以 在 研究 光 变 观测 的 时 间 跨 度 时 估算 小 行星 YORP 效应 的 探测 概率 。 


小 行星 YORP 效应 目前 面临 的 主要 问题 包括 但 不 限于 表面 微 结构 的 量化 和 热 物理 模型 
的 进一步 发 展 ， 未 来 利用 光 变 数据 来 探测 YORP 效应 仍然 是 主要 的 途径 。 随 着 小 行星 的 光 
变数 据 不 断 增 加 ， 通 过 筛选 模型 有 效 地 挑选 可 探测 YORP 效应 的 候选 体 变 得 必要 。 这 有 利 
于 寻找 那些 在 YORP 效应 作用 下 自转 减速 的 小 行星 。 此 外 ，YORP 效应 的 自我 加 强 和 自 


我 限制 过 程 也 是 未 来 需要 研究 的 重点 ， 对 建立 完备 的 YORP 理论 较为 重要 四 
YORP 效应 不 断 发 展 和 完善 的 过 程 中 ， 尚 有 许多 有 意义 的 工作 值得 学 者 们 去 探索 。 
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The Review of the YORP Effect for Asteroids 
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Abstract: When smaller asteroids are heated up by the sunlight, they re-radiate eventually 


the energy away in the thermal waveband, which creates a tiny thrust in turn. Due to irregu- 
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lar shape, thrusts on the different parts of asteroid result in a thermal torque complemented 
by scattered sunlight, which can modify the rotation rates and obliquities of small astreoids. 
This rotational variant has been named as Yarkovsky-O’Keefe-Radzievskii-Paddack (YORP) 
effect. During recent years, the YORP effect has been used to resolve many mysteries in 
asteroid science. In this paper, we review some key results especially in observational field, 
and preview the future goals for this work. 

Firstly, the paper introduces the great significance of YORP research. Beside explaining 
the rotation and size distribution of asteroids, YORP effect also is applied to assessing near 
earth asteroid impact threat precisely, constraining asteroid physical characteristics, explain- 
ing surface material migration and asteroid activity triggering, even act as a important role 
on planetary system evolution. 

Secondly, the paper introduces the theoretical framework of YORP effect, and illustrates 
the principle and application conditions. The torque calculation of the YORP effect and its 
influence on the rotation rate and obliquities of asteroids are reviewed in detail. More similar 
effects, such as BYORP effect, Orbital YORP effect, TYORP effect, are also introduced 
comprehensively. 

Thirdly, the paper presents a progress on the rotational acceleration driven by YORP 
through observation, which is very dependent on the shape inversion model and thermophys- 
ical model of asteroid. So far, only 6 near earth asteroids have been measured and confirmed 
rotational acceleration, and all of them are speeding up. Also, the asteroids which have 
opportunity to measure rotational acceleration in the future are given. 

Fourthly, the paper presents the role of YORP effect in the statistical distribution of 
asteroid rotation state, mainly including: 1) flatting the rotational speed distribution in the 
asteroid belt, and driving that of near earth asteroids as a bimodal distribution; 2) alignment 
of the rotation axes of asteroids, especially evidence in asteroid families. Furthermore, the 
YORP factor (Cy) is used to estimate the degree of YORP effect. 

Finally, the paper discusses some key points should be accounted in target selection 
of measurable YORP effect, also including the influence of surface micro-structure charac- 
teristics, rough surface thermal-infrared beaming and global self-heating of asteroid on the 


quantitative evaluation of YORP effect. 


Key words: asteroid; YORP effect; shape inversion model; thermophysical model 


